Landau mapping and Fermi liquid parameters of the 2D t— J model 
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An exact-diagonalization technique on small clusters is used to study the momentum distribution 
function n(k) of the lightly doped t—J model in two-dimension (2D). We find that n(k) can be 
decomposed into two components with bosonic and fermionic doping dependence. The bosonic 
component originates from the incoherent motion of holes and has no significance for low-energy 
physics. For the fermionic component we explicitly perform the one-to-one Landau mapping between 
all low- lying states of the t—J model and those of a system of spin- 1/2 quasiparticles and extract 
the quasiparticle dispersion and Landau parameters. These results demonstrate that the t—J model 
is a Fermi liquid with a 'small' Fermi surface and a moderate attractive interaction between the 
quasiparticles. 
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Despite considerable experimental efforts, the Fermi- 
surface topology of cuprate superconductors continues to 
pose an intriguing and not really well-understood prob- 
lem. While it seemed to be settled for some time that 
the Fermi surface of these materials is simply the one pre- 
dicted by the local-density approximation (LDA) with an 
only moderate 'correlation narrowing' of the bandwidth, 
recent developments in photoelectron spectroscopy, like 
the discovery of the 'shadow bands' or the tempera- 
ture dependent pseudogap at (it, 0) 0, have challenged 
this point of view. There is moreover the long-standing 
problem that Fermi-liquid-like calculations based on the 
LDA Fermi-surface cannot describe the doping depen- 
dence of either Hall constant or dc-resistivity on the hole 
concentration 6; both quantities consistently suggest a 
carrier density ex S §3], rather than oc (1 — S) as it would 
be for the Fermi surfaces predicted by the LDA. On the 
other hand, taking the shadow bands as true part of the 
Fermi surface, which then would have the topology of 
elliptical hole pockets centered on (tt/2,tt/2), would im- 
mediately lead to complete accord between Fermi-surface 
topology and transport properties in the framework of a 
very simple Fermi- liquid-like picture 0. While such a 
picture thus is quite appealing, clear experimental evi- 
dence for hole pockets has not been found so far H. 
It is the purpose of the present paper to address the prob- 
lem of Fermi-surface topology of cuprates theoretically, 
by studying the 2D t—J model, the simplest strong corre- 
lation model which may give a realistic description of the 
Cu02 plane. To that end we analyze the electron mo- 
mentum distribution (EMD) function n -(k) = (cj C(T CkCT}- 
In both the conventional Fermi liquid, as well as more 
exotic quantum liquids such as the Luttinger liquid, the 
EMD shows singularities which mark the location in k- 
space of the extreme low-energy single-particle-like ex- 
citations. While u(k) (for simplicity we omit the spin 
index unless it is indispensable) can be readily evaluated 
exactly for small clusters of strong correlation models, 



one may not necessarily gain much from such data. The 
reason is that in the photoemission spectra of strongly 
correlated electrons the larger part of spectral weight 
resides in incoherent continua whereas the quasiparticle 
peak near the Fermi energy Ep carries only a small frac- 
tion of the weight. Since rt(k) is the zeroth moment of 
the spectral weight, it can be decomposed as a sum of 
the quasiparticle part, which is responsible for the step 
of magnitude Z at the Fermi surface, and an incoherent 
background: 



n(k) = Z • 6(w Q p - E F ) + /dw J 4 inc (k,w), 



(1) 



where wqp denotes the quasiparticle dispersion and 
Aj nc (k, u>) the incoherent high-energy part of the spectral 
function. In general the latter part has a k-dependence 
of its own which is superimposed over the true Fermi 
surface discontinuities. In principle, the incoherent back- 
ground should be a smooth function of k, so that a suf- 
ficiently good momentum resolution would allow to dis- 
tinguish the step- like variations at Ep; this however is 
out of the question with the relatively coarse k-meshes 
available in cluster diagonalization, so that any simple 
assignment of a Fermi surface is impossible (see Rcf. || 
for a counterexample). On the other hand, numerical 
diagonalization yields the exact EMD for in principle all 
low-energy states without any statistical error, and in the 
following we want to take advantage of this feature. We 
will show that the smoothly varying incoherent compo- 
nent of n(k) is practically identical for all the low-energy 
states with given hole number. This allows to subtract 
off this incoherent component from the EMD for the vari- 
ous low-energy states, and thus make visible the 'coherent 
structures' in n(k), which originate from the true Fermi- 
surface physics. These coherent structures, which give 
the k-space distribution of the quasiparticles rather than 
the bare electrons, then can be used to establish the Lan- 
dau mapping between the exact low-energy eigenstates of 
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the t—J model to those of a suitably chosen quasiparti- 
cle Hamiltonian and, by matching the excitation energies 
of the two systems, we can directly estimate the Landau 
parameters of the t—J model. 

Let us first note that at half-filling the EMD for the 
t—J model is n(k)=l/2, i.e., independent of momentum, 
and that this is similar to a band insulator where we 
would have n(k)=l. Also, the EMD for a state with 
one hole at momentum ko and total spin a® would be 
n CT (k)=l— <5k,-k ' ^cr.ffo m a band insulator, i.e., a con- 
stant with a 'dip' at — ko and do- The true single-hole 
EMD for the t—J model however is quite different as we 
can see in Fig. [|. To emphasize small differences, this 
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FIG. 1. EMD for the single-hole ground state of the 20-site 
cluster (with total momentum ko=(7i"/5, 3tt/5), S z = 1/2, and 
J/t=0.5). The figures show (a) W|(k), (b) nj,(k), and (c) 
ni a c(k), in the entire Brillouin zone. The edge of the square 
centered on a given k-point is proportional to |n CT (k) — N„/N\. 
Positive (negative) values are indicated by black (crossed) 
squares, and the 'calibration' for their magnitude is given in 
the bottom right figure. 

figure shows n (T (k) with its mean value N a /N being sub- 
tracted (where N a denotes the number of cr-spin electrons 
and N is the cluster size). We then find that the dip ex- 
pected at — ko and do does indeed exist but that it is quite 
shallow. Moreover, for both spin directions n CT (k) ac- 
quires an additional, apparently smooth component, and 
the oscillation of this component around the mean value 
is almost independent of a. Another surprising result 
can be obtained by plotting the maximum amplitude of 
the EMD, i.e., A=n(0, 0)— n(n, 7r) as a function of 6 (see 
Fig. ||). We find that up to S as large as 0.2 the relation 
A cx 6 holds with high accuracy, which is on one hand 
an unexpected result if we adopt a free-electron-like pic- 
ture where a change of the electron density should affect 
n(k) predominantly near the noninteracting Fermi sur- 
face (which is far from either (0,0) or (tt, 7r)). We could 
on the other hand understand this result if we were to 
assume that each hole 'brings with it' the same smooth 
contribution as for the single hole (see Fig. pi), so that 



for each hole number n(k) contains this component times 
the number of holes. Such an assumption in turn is very 
natural M if we assume that the doped holes are 
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FIG. 2. Amplitude A=n(0, 0)— n(ir, n) as a function of the 
hole concentration S evaluated for various cluster sizes TV and 
hole numbers Nh, where for odd Nh the average over spin 
directions has been taken. J/t=0A is used. The full line is a 
guide to the eye. 

'spin bags' where the bare hole oscillates rapidly inside 
a region of reduced spin correlation. In this picture the 
smooth component of ra(k) comes from the (incoherent 
and high-energy) oscillation of the hole 'inside' the spin 
bag 0. If this high-energy motion is sufficiently decou- 
pled from much slower drift motion of the entire quasiholc 
(i.e., the hole plus its dressing region), one may expect 
that this high-energy oscillation contributes identically to 
n(k) for each quasihole, whence the scaling of n(k) with 
hole number follows naturally. In 'spin-charge-separation 
language' one might say that the physical quasiparticle 
is a firmly bound state of spinon and holon. The rapid 
oscillation of the bound holon then might correspond to 
the motion of a boson in a localized orbital, and it is easy 
to understand that bosons which populate a localized or- 
bital which is smaller than the cluster size have a mo- 
mentum distribution proportional to the boson (=hole) 
density. Finally, one may expect that the picture pro- 
posed by Lee et al. ||, who modeled the charge degrees 
of freedom by bosons diffusing in a fluctuating magnetic 
field, would also give a diffuse component in n(k) pro- 
portional to the hole number. 

With these ideas in mind we now assume that the smooth 
component in n(k) indeed accurately scales with hole 
number throughout the Brillouin zone; in other words, 
we write the EMD for a single hole: 



r(k) = 



N 



Z<5k,-k • 5a,, 



=(k). 



(2) 



Then, in the first step, we 'close the dip' at — ko by re- 
placing this n(k) value by that for a symmetry-equivalent 
k-point; e.g., in Fig. [l] where k = (7r/5, 3n/5), we replace 
T M(~ 7r /5, — 37r/5) — ► n^(— 37r/5, 7r/5). Subtracting now 
the constant term N&/N, we obtain the incoherent con- 
tribution 7ii, lc (k), which, as discussed above, should cor- 
respond to the oscillation of the bare hole inside the 
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spin-bag-like quasiparticle. For simplicity we average 
ni nc (k) over spin directions and point-group operations 
(see Fig. |l|c). Next, we assume that also at higher dop- 
ing levels each hole gives a two-fold contribution to the 
EMD: i.e., a 'dip', which marks the total momentum of 
the quasiparticle, and the same n inc (k) as for the single 
hole. If this is correct, the EMD for an eigenstate \^ v ) 
with Nh holes (which we assume to be even) should take 
the form 



n(k) 



N-N h 



Z ■n coh (k)+N h -n inc (k). (3) 



Thereby the coherent part, n co h(k), corresponds to the 
EMD of Nh spin-1/2 fermions which form a state with 
the same total momentum, spin, and point-group sym- 
metry as l^&v). It can be readily obtained numerically 
by subtracting + N h • n inc (k) fr om the exact EMD 

and dividing by — Z . To judge the outcome of this pro- 
cedure we now have to know how the Nh quasiholes dis- 
tribute themselves in k-space for given total momentum 
and spin. In other words, we need to guess the quasipar- 
ticle Hamiltonian. The generic form is 



^(l ,,/',;* , • ./, ,S, • S, ;. (4) 



i J a 



We have chosen a next-nearest-neighbor hopping disper- 
sion, i.e., only the hopping integrals t\\ and £20 between 
the second-nearest (i.e., (1,1)) and third-nearest (i.e., 
(2,0)) neighbors are different from zero. Moreover we 
retain nearest-neighbor density-density and exchange in- 
teractions, V10 an d Jio- In view of the many complicated 
ways in which two 'dressed holes' in an antiferromagnet 
can interact with each other [ p^[ , this may be a very 
oversimplified choice; however this will turn out to be 
sufficient to get reasonable results jH]]. The EMD for 
the different eigenstates of this Hamiltonian are easily 
computed and can then be compared to the 'coherent 
components' obtained by the above subtraction proce- 
dure for the t—J model. 

We should note, before discussing the results, that the 
procedure is based on a whole sequence of quite strong 
assumptions: we assume that n(k) really can be mean- 
ingfully decomposed into two components, as in Eq. (^), 
that the smooth component indeed scales accurately with 
the hole number, and consequently is the same for all 
low-energy states of the t—J model, and finally that a 
one-to-one correspondence exists between the states of 
the t— J model and those of the quasihole Hamiltonian. 
The outcome of this procedure, which is shown in Figs. ^ 
and ^, is however all the more surprising. To begin with, 
it is obvious that the calculated 'coherent components' 
n co h(k) are indeed coherent, i.e., in most cases they are 
sharply localized in momentum space and nearly zero al- 
most everywhere in the Brillouin zone. This shows that 
for most k-points our scaling procedure indeed gives a 



good guess for n(k). We would like to stress again that 
we have subtracted one and the same estimated incoher- 
ent contribution from the exact EMD for all the different 
low-energy states in the t—J model to obtain the n co h(k). 
Next, comparing n co h(k) with the EMD for the model 
system, a remarkable correspondence is obvious. Under 
20 different states shown in Figs. || and ^ only the singlet 
state with k=(7r/3, 7r/3) in the 18-site cluster deviates 
significantly. It should be noted that, with the exception 
of the high-symmetry points (0,0) and (tt, it), Figures || 
and |^ comprise the lowest singlet and triplet 
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FIG. 3. Comparison between the 'coherent component' 
ricoh(k) for all exact low-energy states of the 20-site cluster 
t—J model with 2 holes (circles) at J/t=0.5 and the EMD of 
the quasihole Hamiltonian Eq. (^J) with 2 particles (squares). 
Each graph is labeled by the total momentum and spin of 
the respective eigenstate and shows the entire Brillouin zone. 
The diameters of the circles (edges of the squares) centered on 
each k-point are proportional to the respective n(k) (see the 
lower right edge for the 'gauge'). Momenta for the quasihole 
Hamiltonian are —1 times those for the t—J model. 



states for all the allowed momenta in both 18- and 20-site 
clusters; a similarly good correspondence can also be seen 
for the 16-site cluster, which we do not show here. The 
reason for the deviations in the state with k=(-7r/3, 7r/3) 
and 5*=0 are unclear; this state may correspond to a 
charge-density- wave or stripe- like hole arrangement, as 
can be seen for larger values of J Ej . Finally, the corre- 
spondence between the t—J model and the quasiparticle 
system also becomes apparent in a comparison of the ex- 
citation energies in Fig. [|. While we may not expect 
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to obtain perfect agreement with our rather simplified 
choice of the interaction terms in Eq. (|]) , the relative er- 
rors are mostly < 20% and the 'dispersion' of the excita- 
tion energy is reproduced reasonably well. The optimized 
parameter values in Eq. are in=0.255, £20=0. 15, 
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FIG. 4. Same as Fig. |3j but for 2 holes in the 18-site cluster. 
Vw= -0.6375, and J 10 = 0.15. 

In summary, we have shown that the EMD for all the low- 
energy states of the small-cluster t—J model in 2D can 
be accurately decomposed into an incoherent contribu- 
tion, which has a boson-like dependence on hole density 
and is independent of the total momentum and spin of 
the respective state, 
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FIG. 5. Excitation energies Wex for the eigenstates states 
shown in Figs. ^ and ^: triangles are for the t—J model and 
squares are for the quasiparticle Hamiltonian Eq. (^). The 
numbers in the horizontal axis indicate the panels in Figs. |^ 
and ^, counted from top left, right way, to bottom right. 

and a coherent contribution, which corresponds to that 
of spin-1/2 quasiholes. Using the coherent contributions 



we have directly established the Landau mapping be- 
tween the low-energy states of the t—J model and those 
of a 'quasiparticle Hamiltonian' describing the system 
as a Fermi liquid of spin-1/2 particles corresponding to 
the doped holes. The quasiparticle dispersion is dom- 
inated by the next-nearest-neighbor hopping, as would 
be expected for a dispersion dominated by antiferromag- 
netic spin correlation ]l2] ], but also for a 'spinon' in an 
RVB state [^3|. The interaction between quasiparticles 
is dominated by a nearest-neighbor attraction. Our re- 
sults establish, in our opinion conclusively, that the small 
clusters behave like the finite-size equivalents of a Lan- 
dau Fermi-liquid, in that there is a one-to-one mapping 
between the exact eigenstates of the interacting system 
and those of an interacting system of quasiparticles. The 
most obvious extrapolation to the infinite systems is of 
course that the same quasiparticle Hamiltonian also de- 
scribes the thermodynamic limit. There is good reason to 
believe that the cluster diagonalization quite accurately 
describes short-range processes, which are the dominant 
ones in strongly correlated systems; for the paramagnetic 
regime with its very short-range spin correlations it seems 
moreover plausible that longer-range processes are of lit- 
tle importance for either propagation or interaction of 
the holes, so that we believe that the extrapolation of the 
quasiparticle Hamiltonian from the clusters to the infinite 
system is a quite reasonable approximation. Unless the 
interaction between holes (which may be either the inter- 
action intrinsic to the t—J model or the extra Coulomb 
repulsion) drives the system into a charge-density-wave 
state Q (as may be the case in Lai.s75Sr .i25CuO4), 
the system should thus behave like a Fermi liquid with a 
hole-pocket Fermi surface centered on k=(7r/2, 7r/2). 
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